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Abstract 

Measurements of the polarization angle and orientation of cosmological radio sources 
may be used to search for unusual effects in the propagation of light through the uni- 
verse. Recently, Nodland and Ralston have claimed to find evidence for a redshift- 
and direction-dependent rotation effect in existing data. We re-examine these data 
and argue that there is no statistically significant signal present. We are able to place 
stringent limits on hypothetical chiral interactions of photons propagating through 
spacetime. 
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1 Introduction 



The polarization of radiation emitted by distant radio galaxies and quasars offers a way to 
search for chiral effects in the propagation of electromagnetic radiation. Such objects are 
often elongated in one direction, so that one may define a position angle ip which describes 
the orientation of the object on the sky. Synchrotron radiation can lead to a significant 
linear polarization of the source, and the angle x of fhs plane of polarization may also be 
measured fl]]. (The angle of polarization will typically undergo Faraday rotation, but this 
effect can be removed by using the fact that Faraday rotation is proportional to the square of 
the wavelength.) One can therefore study the relative angle x~4' between the position and 
polarization vectors, keeping in mind that this quantity is only defined modulo 180°. It has 
been found 0, Q that % — is not distributed randomly; there is a large peak at x~4' ^ 90°, 
and a smaller enhancement at x ~" ~ 0°. Since many of these sources are at significant 
redshifts, and therefore very far away, testing whether this relationship is maintained for 
distant sources provides constraints on possible chiral effects on the propagation of light 
through the universe, which could rotate x ~ away from the intrinsic value (that which 
would be measured at the source). 

From a field theory point of view, the simplest such chiral effect arises from a Lagrange 
density 

C = -^F.^F'^'' - (PF.^F^-' . (1) 

Here, F^^^, = d^A^, — d^A^ is the field strength for the electromagnetic field A^, F^'^ = 
(l/2)e^'^''°"Fpo- is the dual field strength, and is a pseudoscalar field which need not be 
fundamental (it can be a function of other fields in the theory). We set c = 1 throughout. 
This Lagrangian is the simplest way to couple a neutral pseudoscalar to electromagnetism 
in a parity-invariant way, and often describes the effective coupling of pseudoscalar particles 
(such as pions or axions) to photons. 

Our interest here is in the case where (f) varies only very slowly over extremely large 
distances. In that case an electromagnetic wave travelling through the background field 
will undergo a rotation in its polarization state which depends on the change in 0; such an 
effect arises in a variety of contexts [H-[18|- In the WKB limit where the length scale for 



variations in is much larger than the wavelength of the photon, the polarization angle x 
obeys the simple relation 

Ax = A0 , (2) 
where A indicates the change between source and observer. (Here, and in Eq. (^) below. Ax 
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is measured in radians; elsewhere we measure all angles in degrees; no confusion should arise.) 
This effect is independent of wavelength, and can therefore be distinguished from ordinary 
Faraday rotation. Carroll, Field and Jackiw suggested that observations of polarized 
radio sources provide a stringent test of such an effect, since they afford an opportunity to 
constrain A</) over a large interval in space and time (see also [0, |2( 



The specific model investigated in set d^(j) = — (l/2)p^, where is a 4-vector whose 



expectation value parameterizes violation of Lorentz invariance (as well as CPT ||2T|). It 
was hypothesized that there exists a preferred coordinate frame, close to the background 
Robert son- Walker frame of our universe, in which d^pu = 0. This implies that the predicted 
rotation of the polarization angle for a source at redshift z is given in terms of the timelike 
component po and the spacelike vector p by 

Ax = -^r{po-pcose) , (3) 

where p = \p \ = {^^^PiPjY^"^ ■, d is the angle between p and the direction toward the source, 
and r is the proper spacelike distance traveled. If we take a flat {k = 1) universe as a 
reasonable approximation, we have 

2 



3Hn 



1 - (1 + zr''/^ 



(4) 



where Hq is the Hubble constant today. Regardless of whether or not one is interested in 
tests of Lorentz invariance, Eq. (|^) is a useful parameterization of potentially observable 
chiral effects. 

In it was shown that the radio galaxies at redshift greater than 0.4, with maximum 
polarizations greater than 5%, were strongly clustered around x — 'ip^ 90°, using a sample of 
galaxies and redshifts obtained from the literature 0, ^ Assuming that the timelike 

component po would be significantly larger than the spacelike part p, the limit 

Po < 1.7 X lO-'^^ho GeV (5) 



was obtained, where ho = ifo/100 km/sec/Mpc. Recently, Nodland and Ralston p4[, using 
the same set of data,[| searched for anisotropic effects such as those that would arise from 
a nonzero spacelike part p in Eq. (^ . Surprisingly, they claimed to find a significant signal 

^As mentioned in |2^, there were a handful of transcription errors in the table published in The 
corrected values are: entry 9, coordinates 0106+72; entry 35, coordinates 0459+25; entry 84, ip = 39°; entry 
124, coordinates 1626+27; entry 144, z = 0.054; and entry 153, z = 0.0244. We are grateful to Borge 
Nodland for informing us of these corrections. 
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in the data. Given the fundamental importance of such a result, we have undertaken a re- 
examination of the data, and present our results in this paper. We conclude that the data 
are most consistent with no effect, contrary to |^|. Our disagreement stems primarily from 



the method used to disentangle the 180° ambiguity in the quantity x ~ "4^-, ^"^^ the use of 
randomly generated data for comparison purposes, as will be shown below. 

As this manuscript was being completed we received a preprint by Eisenstein and Bunn 



25 1 , who come to conclusions similar to those expressed in this paper. 



2 Data 



We consider the same set of data as was used in |jTO[, with the corrections noted above. This 
set includes 160 sources, with redshifts as high as 2.012. The first step is to establish the 
existence of a reliable correlation between x if), so that we may test its behavior for 



distant galaxies. Following [^], we have divided the data into distant and nearby objects, 
with the division drawn at z = 0.3; there are 89 sources with z < 0.3, and 71 with z > 0.3. 
In Figures One and Two, we have plotted histograms of x — fo^' these two sets. We have for 
the purposes of these figures defined x — so that it lies between 0° and 180°, and grouped 
the data into bins which are 10° wide. 

In Figure One, representing nearby galaxies, there appears to be evidence for a narrow 
enhancement at x ~ ~ 0°, and a broader peak at x ~ ~ 90°. However, the correlation 
is evidently not very strong. The most likely explanation for this fact is that many of these 
galaxies are members of a different population than those at high redshift, which can be 
observed only if they are of high luminosity, and the lower-luminosity galaxies demonstrate 
a weaker correlation between polarization and position angle. For the more distant galaxies 
shown in Figure Two, there is a very clear peak at x ~ ~ 90°. There is no noticeable peak 
at X ~ ~ 0° ill this sample; again, this may be explained if the galaxies with x " ~ 0° 
are members of a lower-luminosity population. (See |0, ^ for discussion of these correlations 
and their interpretation in terms of models of the sources.) 

It is reasonable to suppose that the galaxies with a higher degree of maximum polarization 
would show any effect more strongly than those that are polarized only weakly. We therefore 
show in Figures Three and Four the same plots as in Figures One and Two, this time limited 
only to those sources with maximum polarization greater than or equal to 5%. These are 
the sources that were analyzed in flOl . 



The peaks found in the previous plots, at 0° and 90° at low redshifts and more dramat- 
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Figure 1: Histogram of number of galaxies vs. x — "0; for galaxies with z < 0.3. 
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Figure 2: Histogram of number of galaxies vs. x — ■0) for galaxies with z > 0.3. 
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Figure 3: Histogram of number of galaxies vs. % — V', for galaxies with 2; < 0.3 and maximum 
polarization > 5%. 
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Figure 4: Histogram of number of galaxies vs. % — 'i/^, for galaxies with 2; > 0.3 and maximum 
polarization > 5%. 
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ically at 90° in the high-redshift sample, are also evident in Figure Three and Figure Four, 
arguably more convincingly. However, any increased correlation is offset to some degree by 
the smaller number of data points. Therefore, in the remainder of this paper we will not 
discard those sources with maximum polarization < 5%; we will thus use precisely the same 



data as were analyzed in |24 



The crux of our disagreement with can be found in Figure Two, the distribution of 



X — fo^^ sources with z > 0.3. As noted in ||I0[, this plot constitutes vivid evidence that the 
polarization angle in these sources is intrinsically perpendicular to the position angle on the 
sky. If the claim of is true, it is necessary to believe that the peak at 90° is an accident. 



and these data are actually drawn from a distribution which is intrinsically centered at 0°, 
with position- and redshift-dependent contributions of order 180°. We will argue in the next 
section that this is not the simplest interpretation of these data. 

3 Constraints on chiral effects 

Searching for a signal in the polarization data is complicated by the fact that x ~ is 
only defined modulo 180°. In testing any specific hypothesis, it is necessary to choose some 
reasonable procedure for resolving this ambiguity. The method chosen by the authors of 



was the following: for any choice of direction for the vector p, define an angle j3 = x— ^^±180° 
which is between 0° and 180° if cos^ > and between —180° and 0° if cos^ < 0, where 6 
(which they called 7) is the angle between p and the direction toward the source. 



It was noted in |2J] that this procedure necessarily introduces correlations between P and 
r cos 6. It would be illegitimate, therefore, to take a statistical correlation between these two 
quantities as itself evidence of a signal in the data. However, if the degree of correlation 
were much higher than that which would be expected if there were no signal in the data, we 
might conclude that there was a measurable effect. 

It is at this point in the analysis that we find two important flaws in the procedure 
followed in EM. First, one must reliably determine the zero point for x~'4'> which would be 



observed in the absence of any chiral effects. In the authors searched for a best flt to 



the data of the form (3 = {l/2)A~^rcos6 + 6, where in the notation of Eq. |^, = p^^. They 
found that the favored value for the zero point was 5 ~ 0°. This seems to be inconsistent 
with the evidence of Figure Two, which exhibits a peak at 90°. The resolution is simply the 
fact that the deflnition of (3, as described above, separates the data into two groups, one 
with —180° < (3 < 0° and one with 0° < f3 < 180°. With this procedure the favored value 
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for 6 will always be near 0°; it arises essentially from taking the average of a group of points 
clustered around 90° and another clustered around —90°. This method of resolving the 180° 
ambiguity is therefore inappropriate for data which lie naturally in the vicinity of 90°. 

Nevertheless, argues that the correlation found is statistically significant, as it was 
only very rarely reproduced in artificially generated sets of data. The procedure for gener- 
ating these sets is the second important flaw that we find. Figure Two provides evidence 
that, regardless of the position of the source on the sky, x — is distributed approximately 
in a Gaussian distribution centered on 90°; a best fit to the Gaussian yields a dispersion of 
a = 33°. Therefore, in searching for position-dependent effects, it is appropriate to compare 



the actual data to data which is generated by drawing from a similar distribution. In 



on the other hand, artificial realizations were generated completely randomly, i.e. from a 
fiat probability distribution for x ~ i^- This has a dramatic effect on the claimed signifi- 
cance of the result. We performed an independent analysis,^ using two different methods of 
generating the artificial data sets: first by drawing from a fiat distribution, and then from a 
Gaussian with the appropriate width. The numbers generated were values of x ~ for the 
positions and redshifts of the 71 sources in the sample with z > 0.3. In 1000 realizations of 
the data drawn from a fiat distribution, in only 7 trials was the significance of the correlation 
greater than that in the actual data; this is comparable to the 6 out of 1000 reported in [p4| , 
and if reliable would be evidence of the existence of a signal. On the other hand, in 1000 
realizations of the data drawn from the appropriate Gaussian distribution, the artificial data 
was more strongly correlated with the hypothesized test function in 911 out of 1000 trials. 
Even if there were no signal at all in the data, we would expect the artificial realizations 
to have a stronger correlation approximately 50% of the time; the fact that our trials had 
better correlations over 90% of the time is due to the fact that the Gaussian slightly under- 
estimates the number of data points near 0°. This result, however, vividly demonstrates our 
main point: the existence of a real enhancement of % — near 90° leads to a spuriously large 



correlation coefficient if one uses the procedure described in p4[. When this enhancement, 
which is consistent with conventional models of the sources, is taken into account, there is 
no sign of an additional effect such as that in Eq. |^. 

There is another way of quantifying our claim that a random distribution centered around 



^ There are two procedures described in |2J]. In the first, the correlation in the artificial data sets was 
calculated using the best-fit direction obtained fi'om the real data. As the hypothesis being tested is that 
there is some direction of anisotropy, not that a specific direction is picked out a priori, it is more appropriate 
to compare correlations in each set of data (real and artificial) calculated using the best-fit direction for that 
set. This is equivalent to the second procedure in p4l, and is the method followed here. 
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90° is a better fit to the data than the correlation proposed in Figure Five is a plot of 
X — "^A as a function of r cos 9, where 9 is defined using the best-fit direction quoted in |^ 
and X ~ is defined to be between 0° and 180°. We may think of this graph as being defined 
on a cylinder, where 0° is to be identified with 180°. With this in mind, we have plotted two 
possible relationships, a solid horizontal line at 90° and a dashed line at (l/2)A7^r cos ^ + 5, 
where we have measured the parameters and 5 from Figure 1(d) of If the relationship 



claimed in ||2^ is correct, the dashed line wrapping around the cylinder should be a better 
fit to the data than the solid horizontal line. This can be measured by calculating 

- 1: (1)^ . («) 

where we take the average error to be a = 33°, although the precise value is irrelevant 
for purposes of comparison. The quantity A,, which represents the difference between the 
predicted and measured value of x — V", is of course subject to the 180° ambiguity; however, 
we can resolve this ambiguity optimistically for each point, by defining —90° < Aj < 90°. 
Using this procedure, we calculate that the best fit proposed in yields = 161, while 
the hypothesis of no effect yields = 69. Thus, the horizontal solid line in Figure Five is a 
much better fit than the diagonal dashed lines. 

Given that there is ample evidence that the intrinsic zero point is centered on x — ^ = 90°, 
we may ask how good a limit we can place on an effect such as that in Eq. (^. One 
approach to this problem is to define x ~" to be between 0° and 180°, and to assume that 
the deviation from the intrinsic value is given by Ax = X ~ '4^ ~ 90°. It is then possible 
to do a straightforward least-squares fit to Eq. (^, with the four components of as free 
parameters. Using the data at redshifts z > 0.3, the best-fit parameters obtained in this 
way are 

Po = (0.59 ± 0.80)iJo 

= (1.25 ± 1.71) X 10-^2/^0 GeV , 

\p\ = (1.13 ± 1.40)iJo 

= (2.41 ± 2.99) X lO^^^ho GeV . 

(This procedure yields separate values for each of the three spacelike components pf, since 
each value is consistent with no preferred direction, it is more appropriate to quote the limit 
on the magnitude \p\.) These values are consistent with p^ = 0, and similar to the limit on 
Po from quoted in Eq. (|^). 
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Figure 5: The difference between polarization and position angles as a function of rcosO 
for the best-fit direction of anisotropy proposed in |Q. Angles of 180° are to be identified 
with 0°; the data thus live on a cylinder. The solid line represents the predicted relationship 
in the absence of any signal, while the diagonal dashed line wrapping around the cylinder 



represents the model suggested in [24 
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4 Discussion 



After analyzing the data in a variety of ways, we are able to conclude with confidence 
that there is no evidence for a chiral effect on the propagation of photons from distant 
radio sources. Despite this negative result, there are still good reasons to further pursue 
observations such as those examined in this paper. 

In Figures Six and Seven we have plotted the position of the sources in the sky, indicated 
by symbols related to the deviation of x — from 90°. Figure Six includes all of the galaxies, 
while Figure Seven is limited to the distant sources with z > 0.3. The squares represent 
sources with x ~ < 90°, while the x's are sources with x ~ 4' > 90°. The size of the 
symbol is related linearly to the deviation from 90°, although for clarity there is an offset 
so that points with x ~ very close to 90° still have a nonzero size. One conclusion to 
be drawn immediately from these graphs is that there is a need for additional data to be 
collected in the southern celestial hemisphere, especially at high redshifts. In the future, 
observations of polarization of the cosmic microwave background may be the best source of 



data for constraining phenomena such as these [|T^, p^ . 

In characterizing the limits one can place on chiral effects, for convenience we hypothe- 
sized a fixed four-vector which would represent a violation of Lorentz invariance. If an 
effect were to be found, however, it is by no means necessary that such a profound conclusion 
would have to be drawn. A more plausible hypothesis would be that of a very slowly- varying 
scalar field with a coupling as in the application of the data discussed in this paper 
to this possibility was examined in [|12|. Such a field could arise as an ultralight axion, with 
mass of order the Hubble constant today (or less). Interestingly, such axions may appear 
naturally in the strongly coupled limit of heterotic string theory ^ . Another possibility 
is the detection of axion-like cosmic strings; in the vicinity of such a string, the polarization 
angle of two light rays passing on either side will undergo rotations in opposite directions 0. 
Although there is no obvious sign of such a signal in Figures Six and Seven, the importance 
of such a finding encourages us to continue the search. 
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Figure 6: Positions of radio sources on the sky, including all redshifts. The symbols indicate 
deviations from x~'^ = 90° j squares are sources with x — V' < 90° i x 's are sources with 
X — -0 > 90°. The size of the symbol indicates the amount of deviation from 90°. 
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Figure 7: Same as Figure Six, including only galaxies with z > 0.3. 



15 



5 Acknowledgments 

We would like to thank liana Harms, Jack Hughes, Arthur Kosowsky, and Uros Seljak 
for helpful conversations, and Borge Nodland for assistance with the data. This work was 
supported in part by the National Science Foundation under grant PHY/94-07195 and NASA 
under grant NAGW-931. 

References 

[1] D.J. Saikia and C.J. Salter, Annu. Rev. Astron. Astrophys. 26, 93 (1986). 

[2] J.N. Clarke, P.P Kronberg, and M. Simard-Normandin, Mon. Not. R. Astr. Soc, 190, 
205 (1980). 

[3] P. Haves and R.G. Conway, Mon. Not. R. Astr. Soc, 173, 53p (1975). 

[4] P. Sikivie, Phys. Lett. B137, 353 (1984). 

[5] L. Maiani, R. Petronzio, and E. Zavattini, Phys. Lett. B175, 359 (1986). 

[6] C. Wolf, Phys. Lett. A132, 151 (1988). 

[7] S.G Naculich, Nucl. Phys. B296, 837 (1988). 

[8] A. Manohar, Phys. Lett. B206, 276 (1988); (E) B209, 543 (1988). 

[9] J.A. Harvey and S.G. Naculich, Phys. Lett. B217, 231 (1989). 

[10] S.M. Carroll, G.B. Field, and R. Jackiw, Phys. Rev. D 41, 1231 (1990). 

[11] C. Wolf, Phys. Lett. A145, 413 (1990). 

[12] S.M. Carroll and G.B. Field, Phys. Rev. D 43, 3789 (1991). 

[13] D. Harari and P. Sikivie, Phys. Lett. B289, 67 (1992). 

[14] B.A. Campbell, N. Kaloper, R. Madden, and K.A. Olive, Nucl. Phys. B399, 137 
(1993); |hep-th/9301129| . 

[15] S. Mohanty and S.N. Nayak, Phys. Rev. Lett. 70, 4038 (1993); (E) 71, 1117 (1993); 
(E) 76, 2825 (1996); |hep-ph/93033T0i 

16 



[16] E.D. Carlson and W.D. Garretson, Phys. Lett. B336, 431 (1994). 

[17] L. Masperi and S. Savaglio, Astropart. Phys. 3, 209 (1995); |astro-ph/9411047 . 

[18] M.P. Haugan and T.F. Kauffmann, Phys. Rev. D 52, 3168 (1995); kr-qc/9504032 . 

[19] A. Cimatti, S. di Serego Alighieri, G.B. Field, and R.A.E. Fosbury, Astrophys. J. 422, 
562 (1994). 

[20] M. Goldhaber and V. Trimble, J. Astrophys. Astr. 17, 17 (1996). 
[21] S. Coleman and S.L. Glashow, HUTP-97/A008; |hep-pli/ 9703240 . 
[22] G. Burbidge and A.H. Crowne, Ap. J. SuppL, 40, 583 (1979). 

[23] H. Spinrad, S. Djorgovski, J. Marr, and L. Aguilar, Pub. Astron. Soc. Pacific, 97, 932 
(1985). 

[24] B. Nodland and J.P. Ralston, Phys. Rev. Lett. 78, 3043 (1997); |astro-pli/970IT^ . 

[25] D.J. Eisenstein and E.F. Bunn, ^stro-ph/ 9 704247 . 

[26] U. Seljak and A. Kosowsky, private communications. 

[27] T. Banks and M. Dine, Nucl. Phys. B479, 173 (1996); |hep-th/9605T36 . 

[28] T. Banks and M. Dine, hep-th/9608T97 . 



17 



